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ABSTRACT 

We report the results of a search for an emission line from radiatively decay- 
ing dark matter in the Chandra X-ray Observatory spectrum of the ultra-faint 
dwarf spheroidal galaxy Willman 1. 99% confidence line flux upper limits over 
the 0.4-7 keV Chandra bandpass are derived and mapped to an allowed region in 
the sterile neutrino mass-mixing angle plane that is consistent with recent con- 
straints from Suzaku X-ray Observatory and Chandra observations of the Ursa 
Minor and Draco dwarf spheroidals. A significant excess to the continuum, de- 
tected by fitting the particle-background-subtracted source spectrum, indicates 
the presence of a narrow emission feature with energy 2.51 ± 0.07(0.11) keV and 
flux [3.53 ±1.95(2.77)] x 10'^ photons cm'^ s'^ at 68% (90%) confidence. Inter- 
preting this as an emission line from sterile neutrino radiative decay, we derive 
the corresponding allowed range of sterile neutrino mass and mixing angle using 
two approaches. The flrst assumes that dark matter is solely composed of sterile 
neutrinos, and the second relaxes that requirement. The feature is consistent 
with the sterile neutrino mass of 5.0 ± 0.2 keV and a mixing angle in a nar- 
row range for which neutrino oscillations can produce all of the dark matter and 
for which sterile neutrino emission from the cooling neutron stars can explain 
pulsar kicks, thus bolstering both the statistical and physical signiflcance of our 
measurement. 
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Introduction 



1.1. Context 



None of the Standard Model particles can account for the dark matter that makes up 
most of the mass in the universe. We report the latest results of a search for dark matter in 
the form of relic sterile neutrinos (see Kusenko 2009 for an up-to-date review). The motiva- 
tion for considering the sterile neutrino, one of a number of feasible dark-matter candidates, 
is two- fold. First, the discovery of ordinary neutrino masses is most easily accommodated by 
means of the so-called seesaw mechanism, which calls for some new gauge-singlet fermions. 
If all of these fermions have very large Majorana masses, there are no additional degrees 
of freedom (particles) at the low energy scale. However, if one of these mass parameters 
lies in the 1-30 keV r ange, the corresponding st erile neutrino can be the long sought-after 



dark matter particle ( jPodelson fc Widrowlll994l ). Second, the physics of supernovae, that 



often assists in ruling out hypothetical low-mass particles, provides some intriguing clues 
in favor of the existence of sterile neutrinos with the same parameters that are required to 
explain dark matter. Such particles would be anisotropically emitted from a cooling newly 
born neutron star, inducing a suf ficient recoil mornentum in the neutron star to explain 
the observed velocities of pulsars ( iKusenko fc Segre I1997I : iFuUer et al.lboOsL X-ray obser - 



vations offer the best probe of this well-motivated dark-matter candidate (jKusenko 2009). 



(Dodelson & Widrow 


199^ 


1), resonant oscillations ( 


Shi & Fuller 


1999) and via various other 


channels ( 


Kusenko 


200d: J 


^haposhnikov & Tkachev 


2006|; 


Petraki & Kusenko 


2008 


). The ki- 



netic properties of dark matter relevant for structure formation on small scales (i.e., how 
"warm" or "c old" the dark matter is) depend on the produc tion scenario, as well as the 



particle mass ( iKusenkd l2006l : iPetrakil l2008l : iBoyanovskyI l2008bl ) . 



The generi c prediction is that relic sterile neutrinos must decay into a lighter neutrino 
and a photon (jPal fc Wolfensteiru Il982l ). Since this two-body decay of a 1-30 keV mass 
particle produces a narrow line with energy Ej = mstC^/2, X-ray astronomy provides a 
unique opportunity to discover these relic particles. Archival X-ray data has been used to 
set limits on relic sterile neutrinos (see the review in Kusenko 2009). In addition, a dedicated 
sea rch for sterile neutrinos using the Suzaku X-Tay telescope has recently been conducted 
by iLoewenstein. Kusenko. fc Biermanru (120091 ) (hereafter, LKB). Here we will present new 
results from the search using the Chandra X-ray Observatory. Our analysis of data on 
the Willman 1 dwarf spheroidal galaxy yields evidence of a 2.5 keV line, consistent with the 
decay of a relic sterile neutrino with mass 5 keV; and, with an inferred line strength consistent 
with that expected if all of the dark matter is composed of sterile neutrinos produced by 
oscillations. 
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1.2. Target selection 



Given tlie uncertainty in sterile neutrino properties, and the importance of dark matter 
discovery, one ought to explore the X-ray band with existing instrument s using well-chosen 
targets. While the dec ay timescale greatly exceeds the Hubble time (jPal fc Wolfenstein 
19821 : iBarger et al.lll995l ). the emissivity may reach detectable levels in regions of large dark 
matter s urface density, such as g alaxy clusters, M31, and Local Group dwarf spheroidal 
galaxies (lAbazajian et al.l l2001al Jbl: iDolgov fc HansenI l2002l : iBoyarsky et al.l l2006l ) . The last 
offer, arguably, the most promising opportunity because of the proximity and high dark mat- 
ter density in these systems, and the absence of additional competing internal X-ray sources. 
Suzaku Observing Cycle 2 observing time was awarded time to study the best available 
targets at that time: the Ursa Minor and Draco dwarf spheroidals (LKB). Subsequently, a 
new sub-population of faint Milky Way satellites discover ed in the Sloan Digital Sky Survey 



(SDSS) were found to have very high mass-to-light ratios (IStrigari et al.ll2008al : I Walker et al. 



20091 ) . These rival or surpass Ursa Minor and Draco as compelling targets for sterile neutrino 
searches. We were awarded 100 ks of observing time in Chandra cycle 10 to investigate one of 
these - namely Willman 1. These data represent the first X-ray observations of this extreme 
system; and, we report on our findings here. 



1.3. Previous Results on Ursa Minor 

The basic equations relating the galaxy parameters (distance lOOrfioo kpc, dark matter 
mass in projection Mpro = lO^My M©), X-ray observables (line energy E^, line flux Fiine), 
and sterile neutrino parameters (mass rrist, mixing angle 6, fraction of dark matter in sterile 
neutrinos fst), are the following (Loewenstein et al. 2009, Kusenko 2009, and references 
therein) : 

m^t = 2E^, (1) 



(2) 



and 



Fn,e = 4.69 X 10-6 r_27/st Mr rfroo (^^j (3) 
= 5.15 X 10-1° (^) /st Mr d^ol photons cm-^ (4) 



The abundance of sterile neutrinos produced through non-resonant oscillations may be cal- 
culated for given values of nisi and 6, assuming the standard early thermal history of the 
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universe. This relation determines fst as a function of m at and 9, assuming that all sterile 
neutrinos are produced in this way as first proposed by iDodelson fc Widrowl ( 19941) . and 
lower limits on f^t if there is some non-negligible lepton asymmetry ( Shi &: Fuller Il999 ) or if 
some other prod uction channels contributed more than the oscillations to the sterile neutrino 
relic abundance (lKusenkdl2009l ). 

LKB did not detect a sterile neutrino emission line in Ursa Minor Suzaku spectra. The 
line fiux upper limit was used to derive an excluded region in the rrist — 9 plane for the case 
where all of the dark matter is composed of sterile neutrinos. A more general constraint 
was inferred from the condition that the line fiux upper limit not be overproduced through 
oscillations at the sterile neutrino relic density corresponding to each pair {mst,9). We 
derived line fiux upper limits using a maximum-likelihood approach that, along with the 
lack of intrinsic X-ray emission, enabled us to minimize systematics, and account for those 
that remained. The resulting constraints on sterile neutrinos are illustrated in Figure 9 
of LKB. Although the constraints from this single Suzaku observation are comparable to 
previous results over the entire 1-20 keV mass, the allowed range does not rule out sterile 
neutrinos as a viable dark matter candidate. We therefore expanded our search to the 
Willman 1 system - conducting the first X-ray observation of one of the new SDSS dwarf 
spheroidal galaxies. The capabilities of the Chandra X-Ray Observatory are well-suited to 
the compact size of Willman 1; and, these data provide a good complement to the broader 
beam observations of the more extended classical Milky Way dwarf spheroidals (Section 4.2). 



2. Data Analysis 

2.1. Observation and Data Processing 

Willman 1 ( Chandra ObsID 10534, SeqNum 600727) was observed by Chandra on Jan- 
uary 27-28 2009 for 102745.5 seconds. The data was taken in TIMED VFAINT mode, and 
was pipeline-processed on 2009-2-10 using ascdsver 7.6.11.10. We use CIAO version 4.1.2 
(CALDB 4.1.4), and XSPEC version 12.5 in our analysis. Chips 0, 1, 2, 3, 6, 7 were on 
during the observation. Willman 1 was positioned at the ACIS-I aimpoint, and we generally 
restrict our analysis to the I-array (chips 0-3) data. 

We analyze the standard level 2 event file, and also reprocess the level 1 event file in 
order to implement the VFAINT data-mode quiescent background reduction techniquejl] For 
the latter, likely background events are fiagged based on the pulse heights in border pixels of 



http://cxc.harvard.edu/ciao/threads/aciscleanvf 
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the 5x5 pixel event islands; default split threshold and trail parameters are adopted. These 
events, as well as those corresponding to bad grades, are removed as the new level 2 event file 
is generated. No flares are found upon examination of the light curves; and, time intervals 
where the count rate deviates by more than 10% of the mean are excised. The final exposure 
times are 99400 and 97100 seconds for the cleaned pipeline processed and reprocessed event 
files, respectively. 

2.2. The Background 

We consider off-axis (S2, i.e. chip 6), blank sky, and stowed data sets as candidates for 
constructing background spectra. The last corresponds to events collected with the ACTS 
stowed while the HRC-I is in the focal plane; and, has been demonstrated to be consistent 
with the dark moonj^ That is, it can provide an estimate of the pure particle background 
(PB). In contrast, the blank sky and S2 data include the cosmic (CXB) and galactic (GXB) 
X-ray backgrounds (which vary over the sky, especially the GXB). The blank sky and stowed 
event files are reprojected using the Willman 1 aspect solution. The quiescent PB reduction 
procedure may also be applied to these background event files. If dark matter is composed of 
sterile neutrinos, their presence in the Milky Way results in line emission in the S2 and blank 
sky data that will reduce the signal if subtracted out as background. Additionally, the S2 chip 
is still within the solid angle subtended by the Willman 1 dark halo (see below). Therefore, 
our constraints are derived from the PB-subtracted spectra from which no astrophysical 
sources of extended X-ray emission have been eliminated. 

2.3. Spectral Analysis Setup 

Source and background spectra are extracted from identical 5' circular apertures cen- 
tered on the Willman 1 position or its equivalent in detector coordinates, except for the S2 
chip background spectrum that is extracted from a 3' circle optimized in a way that accounts 
for chip uniformity and point sources near the edge of the chip. The S2 chip center is 15' 
(165 pc) from the position of Willman 1, and lies outside the region where stellar velocities 
are measured. The allowed dark matter profile thus spans a large range: the average surface 
mass density in this field may be as large as ~ 50%, or as small as ~ 10%, that in the central 
I-array extraction region. Regions corresponding to point sources detected in the Willman 
1 image via application of the CIAO wavelet point source detection algorithm are excluded. 



^http://cxc. harvard.edu/contrib/maxim/stowed/ 
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Background spectra are re-scaled based on the relative count rates in the 9-12 keV band 
where virtually all events are due to the PB. The PB has been found to vary significantly in 
normalization, but not in shape, so that the rescaled PB spectra can be used to construct 
an appropriate PB-subtracted source spectrumjfl 



2.4. Spectral Analysis 

In our analysis of the Suzaku data on the Ursa Minor dwarf spheroidal galaxy, we based 
our constraints on the total (unsubtracted) spectrum. Because the Chandra PB is higher, 
but more precisely determined, we adopt a more conventional approach here and consider 
background-subtracted spectra, grouped to a minimum of 15 counts per bin, utilizing 
statistics to derive best fits and parameter uncertainties. 

The cleanest PB subtraction is realized utilizing the source and stowed background spec- 
tra with application of the VFAINT data-mode quiescent background reduction technique, 
although using the corresponding pipeline-processed datasets yielded similar results (with 
slightly larger uncertainties). The source spectrum is fitted with a model consisting of a 
solar abundance thermal plasma (to model the GXB) and power-law (to model the CXB) 
absorbed by the Galactic column density. In this case the GXB component is not required, 
although the upper limit on the fiux is consistent with the larger-scale average GXB surface 
brightness measured in the ROSAT All-Sky Survey. The best-fit {x^/v = 1.09) power-law 
photon index is 1.8 (1.4-2.7 at 90% confidence), and the normalization consistent with the 
fraction of the CXB unresolved by our observation (~ 25%; Mushotzky et al. 2000). Figures 
la and lb show the total, (stowed) background, and background-subtracted spectra for cases 
utilizing pipeline-processed, and reprocessed (reduced-background) event lists. The best-fit 
power-law model, and fit residuals, are also shown. We show the full 0.3-12 keV energy 
band, although the fits are conducted over 0.4-7 keV. There are PB-related residuals in the 
background-subtracted spectrum above 7 keV where the Chandra sensitivity experiences a 
steep decline. Some of these persistent residuals may be attributed to the fact that the 
time-dependent gain correction cannot be re-applied to the background data, as was done 
with the source data. 

The blank-sky background (SB) is significantly brighter than the source spectrum below 
1 keV. The total PB-subtracted 0.4-7 keV surface brightness is ~ 1.7x higher in the SB than 
in Willman 1. Evidently the GXB is brighter in the former and, as a result, the SB-subtracted 
spectrum ought to be utilized with care (Figure 2). Also, as previously mentioned, if there is 



•^http: / /cxc. harvard.edu/contrib/maxim/bg/ 
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Fig. 1. — (a) top: Total (red), PB (green), and source {i.e. PB-subtracted; black) spectra. 
The histogram is the best-fit GXB+CXB model to the source spectrum, (b) bottom: 
Same as top with the quiescent-particle-background reduced in the total and PB spectra as 
explained in the text. 
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dark matter decay in the X-ray band there will be a contribution from the Milky Way halo 
to the SB that is absent in the PB. 



3. Limits on Sterile Neutrino Parameters from the Chandra Spectrum of 

Willman 1 

In deriving upper limits, we divide the spectrum into low (< 1.1 keV) and high (> 1.1 
keV) energy segments, since any GXB emission will be negligible in the latter. The high 
energy segment best-fit {x^/u = 1.11) power-law slope is 1.77 (1.22-2.37 at 90%, confidence). 
At low energies, power-law jv = 0.95, slope = 2.27^J;3g) and thermal apec {x^ jv = 0.92, 
kT = 0.48 ± 0.15 keV) models provide acceptable fits. We conduct Monte Carlo numerical 
experiments (LKB, and references therein) where we find that, in fits to simulated continua 
spectra with models that include an emission line, the fit is improved by Ax^ = 9-2 with ~ 1% 
frequency - closely approximating that for a distribution with two degrees of freedom. We 
thus adopt Ax"^ = 9.2 fiux limits, derived from adding an unresolved Gaussian component 
in 10 eV intervals over the 0.4-7 keV bandpass to simulated spectra of the best-fit line-free 
model. The errors are averages over 100 simulations that are initiated using parameters from 
fits to a seed simulation (and not the best-fit parameters to the actual data), so that the 
uncertainties in the parameter values are taken into account. Line fiuxes are allowed to go 
negative (as is the case with the simulations). The upper and lower line fiux limits, in 200 
eV bins, based on the PB-subtracted fits are shown in Figure 3. The latter are negative over 
the entire 0.4-7 keV bandpass. The constraints based on the SB-subtracted fits are similar 
in magnitude above 2 keV. 

The mapping of the Willman 1 constraints onto the sterile-neutrino-mass-mixing-angle 
plane (equation 4) is shown in Figure 4a. As discussed above (Section 1.3; see also Section 
3 in LKB and Section 5 below) there are two kinds of limits, both of which are shown. 
The first assumes only that the thermal history of the universe is given by standard Big 
Bang cosmology. Sterile neutrinos are ruled out for the region in parameter space where 
the flux from sterile neutrinos prod uced at minimal abundance via the DW mechanism 



(lAsaka. Laine. fc Shaposhnikovl 120071 ) exceeds our derived limits. The second limit is that 
which overpredicts the derived limits if sterile neutrinos account for 100% of the dark matter 
(dark-matter fraction of sterile neutrinos f^t = 1), regardless of how they are produced. In 
Figure 4b, the former limits, which are the most general, are compared with those we derived 
in Ursa Minor with Suzaku. 



The smaller Chandra beam and more compact Willman 1 dark matter distribution com- 
pensate somewhat for the higher sensitivity of the Suzaku measurement to yield constraints 
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Fig. 2. — Total (red), SB (green), and SB-subtracted (black) spectra. The histogram is the 
best-fit GXB+CXB model to the PB-subtracted spectrum of Figure lb. The SB results in 
an oversubtraction below 2 keV. Due to the contribution of the Milky Way halo the SB is 
expected to contain the dark-matter signal as well; and, therefore subtraction of the SB can 
simultaneously lower both the background and the signal. The oversubtraction below 2 keV 
and the weakening of the signal due to subtraction of the Milky Way contribution are the 
two reasons why we do not use the SB-subtracted spectrum for data analysis. 
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Fig. 3. — Ax^ = 9.2 line flux upper (solid histogram) and lower (broken histogram) limits 
in photons cm~^ s~^. 
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at sterile neutrino masses < 14 keV (corresponding to photon energies in the spectra < 7 
keV) of similar magnitude to those derived for Ursa Minor in LKB. This comparison utilizes 
a projected mass enclosed within the relevant radius of 55 pc (5' at 38 kpc) Mpj-o = 2 x 10^ 
Mq - corresponding to a m ass surface density of S == 210 MqPc"^ - as determined by 



projecting the best-fit NFW (INavarro. Frenk. fc White! 119971 ) mass model in IStrigari et al. 



(l2008bl ). The 90% confidence contours for the NFW parameters allow values of Mpro that 
are 4 x lower, and 2 — 3 x higher. In addition, we estimate that the Milky Way dark matter 
halo and the extragalactic dark matter contribute an additional 65 — 125 Mopc^^ (Xue at 
al. 2008, and references therein). While maintaining a fiducial mass of 2 X 10^ Mq, we show 
the effect on the general constraints of varying the surface density, and hence the effective 
mass in the FOV, by factor of two in either direction in Figure 4b. 



4. A Candidate Line at 2.5 keV 



In addition to the limits derived above, we present evidence of a spectral feature con- 
sistent with a 2.5 keV line from decay of a 5 keV relic sterile neutrino. 



4.1. Line Flux and Statistical Significance 



The largest residual excess emission to the PB-subtracted Willman 1 spectrum with 
respect to the best-fit continuum model occurs at ~ 2.5 keV. When we add an additional 
Gaussian component, we find that improves by 5 (Figures 5 and 6). The line-width 
is consistent with an unresolved feature broadened only by the instrumental line spread 
function, as expected considering the small Willman 1 velocity dispersion. For power-law- 
plus-Gaussian fits to the 2.2-5 keV spectrum, the best-fits and uncertainties on the line 
energy and flux are 2.51 ± 0.07(0.11) keV and [3.53 ± 1.95(2.77)] x 10"^ photons cm'^ s-\ 
where the errors are 68% (90%) confidence limits for two degrees of freedom corresponding 
to Ax^ = 2.3 (4.61) (Section 3; Yaqoob 1998). The corresponding confidence contours are 
shown in Figure 7, as is the Ax^ = 9.21 (99% confidence) contour. The 2.2 keV limit was 
chosen to avoid the energy associated of the peak of the brightest of the PB Au fluorescence 
lines, and with the deep Chandra mirror iridium edge (Chartas et al. 2000; there is also a 
shallow edge at ~ 2.55 keV that does not impact our results). The issue of these instrumental 
features is specifically addressed in Section 4.2.4. However, extending the range to the full 
1.1-7 keV high energy segment bandpass yields a consistent measurement (-Fune = [2.85 ± 
1.85(2.57)] X 10~^ photons cm~^ s~^) as do fits over other subsets of this bandpass. In the 
2.4-2.6 keV energy interval, the total, PB, and source counts are 696 ± 26.4, 601 ± 13.6, and 
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Fig. 4. — (a) top: The sterile neutrino parameter space to the right of the sohd curve is 
excluded by the Chandra observation of Willman 1 assuming only that DW production by 
neutrino oscillations takes place. The sterile neutrino parameter space to the right of the 
broken curve is excluded if 100% of dark matter is composed of sterile neutrinos produced 
by some (unspecified) mechanism . The mass of 2 x 10^ Mf;^ is d e rived fr om the best-fit NFW 
(INavarro. Frenk. &: White! Il997| ) mass model in IStrigari et al.l ( l2008bl ). (b) bottom: The 
solid curve from (a) is reproduced and compared to constraints assuming projected masses 
of 1 and 4 x 10^ Mq (broken curves; see text), and to the constraints derived from Suzaku 
observation of Ursa Minor (dotted line). 
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95 ±30, respectively. In the composite models the emission hne component accounts for 90% 
of the flux. 20 counts in this interval are expected based on the best-fit power-law model 
without a line. Figure 8 shows the unfolded spectrum - the ratio of the data values to the 
product of the model multiplied by the response, all multiplied by the model values - over 
the 1.5-5 keV spectral region. 

The feature is not positively detected in the global blank sky spectrum (Abazajian et 
al. 2007, and references therein). The Ax^ = 4.61 upper limit based on the reprojected 
estimated SB of Willman 1 is 2.1 x 10^^ photons cm^^ s^^. When we utilize the SB, 
or appropriately-scaled S2, spectrum as background the fits are not significantly improved 
in models with an additional emission line, although the upper limits are consistent with 
the hne flux based on the PB-subtracted spectrum. As the SB (Figure 2) and S2 spectra 
overestimate the Willman 1 background, and will include a contribution due to the Milky 
Way dark matter halo (plus, possibly significant, additional emission from the outer Willman 
1 halo in the latter), the expectation that the PB-subtracted spectrum is most sensitive for 
these purposes is confirmed. 

Notably, the fiux is consistent with the narrow region in the mass-mixing-angle plane 
where the DW neutrino oscillation mechanism produces all of the dark matter (and where 
sterile neutrino emission from neutron stars may explain pulsar kicks - sec next section). 
Figure 9 shows the predicted flux at each energy assuming DW production of sterile neu- 
trinos of rrist = 2Ey, and fst = 1. We compare this with the Ax^ = 4.61 upper limits derived 
from Monte Carlo simulations (see Section 3), and our measurement. Even accounting for 
uncertainties in the DW production rate, and allowing for a factor of four uncertainty in 
the mass, a fluctuation at this level maps into this narrow region only over ~ 1 — 3.5 keV. 
At lower energies the predicted line flux is below our sensitivity, while at higher energies 
a much stronger feature would be expected. That is, the likelihood of a marginal feature 
compatible with the simplest sterile neutrino hypothesis is much less than the likelihood at 
some random X-ray energy. 



4.2. Discussion 

4-2.1. Comparison with Limits from Classical Milky Way Dwarf Spheroidals 

We re-examine the Suzaku spectra of the Ursa Minor and Draco "classical" Milky Way 
dwarf spheroidals to check for consistency with the line flux limits of the 2.5 kcV candidate 
sterile neutrino radiative decay emission line in Willman 1. The projected mass estimate 
adopted in LKB within the Suzkau spectral source region for Ursa Minor is ~ 30 x our 



- 14 - 



> 



T3 



2.4 2.6 2.8 

Energy (keV) 

Fig. 5. — (a) Data and best-fit power-law model to the 2.2-5 keV spectrum, in the 2.2-3 keV 
energy region (top), and (b) contributions to over 2.2-5 keV(bottom). 
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Fig. 6. — Same as Figure 5 for the best-fit power-law-plus-Gaussian model. 
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Unfolded 1.5-5 keV spectrum; see text for details. 
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Fig. 9. — The histogram shows the Ax^ = 4.61 upper hmits for a hne-free spectrum. The 
sohd and broken Unes show the expected hne flux from sterile neutrino radiative decay in 
Willman 1, assuming fst = 1 and taking into account the hadronic uncertainties in produc- 
tion (Asaka et al. 2007), which affect the relation between the mass and the mixing angle 
corresponding to a given value of fst- The dotted curves allows for a factor of two uncertainty 
in mass in either direction (My = 0.1 — 0.4), in addition to the spread in possible production 
rate. 
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fiducial value for tlie mass witliin tlie (smaller) Willn aan 1 Cfeandra extrac tion region. Draco 



has a similar projected mass to that of Ursa Minor ( IStrigari et al.ll2007a( ). 



Factoring in these mass ratios, and their respective distances, the 90% limits of the 
Willman 1 line flux translates to a predicted flux of (3.2 ± 2.5) x 10~^ photons cm~^ s~^ in 
Ursa Minor. In re-analyzing these spectra, we focus on the 1.3-5.2 keV energy region and, 
as in LKB, fit the unsubtracted and unbinned spectrum to a model that includes power-law 
and Gaussian non- X-ray background (PB) components (see LKB for details on how the 
spectra are segmented, and other details). An additional line is not required based upon 
simultaneous re-analysis of the Suzaku XISl (frontside-illuminated) and XIS03 (co-added 
backside illuminated) spectra - in fits where the line flux is tied in the two spectra. However, 
fits to the XISl spectra, when analyzed independently, are improved by the addition of such 
a line (AC=7, where C is the Cash-statistic; Cash 1979). The best fit line energy is 2.41 
keV with 90% confidence limits 2.35-2.57 keV (and a secondary local minimum at ~ 2.5 
keV); the line flux is (2.5 ±1.9) x 10^^ photons cm"^ s^^ (Figure 10a) - consistent with 
the expectations. However, the line is not detected in the (less sensitive) XIS03 Ursa Minor 
spectrum: the upper limit is 1.8 x 10~^ photons cm~^ s~^ (Figure 10b). The Suzaku spectra 
of Draco, reduced and analyzed in an identical manner, do not require a 2.5 keV line. The 
upper limit from simultaneous fits is ~ 2.5 x 10~^ photons cm~^ s~^, with similar limits 
from the individual detectors. Because of its larger distance, the expected flux from Draco 
is ~ 25% smaller than for Ursa Minor. 



Riemer-Sorensen fc HansenI ( 120091 ) analyzed archival data (ObsIDs 9589 and 9776) of 
the Draco dwarf spheroidal galaxy. We have re-analyzed these data, essentially replicating 
their reduction procedures, but deriving line flux upper limits from the spectrum extracted 
from a 4' circular aperture by application of the procedure used to derive limits for Will- 
man 1 (Section 3). There are significant differences in the two datasets. The Willman 1 
observation is about three times deeper; and was centered on the I-array instead of the S3 
chip. The latter has ~ 30% higher effective area at 2.5 keVfl but also a higher background. 
Other factors that contribute to a lower background in the Willman 1 observation are the 
lower point source detection threshold (about 5 times as many point sources are detected), 
the higher concentration of Willman 1 that allows us to use a a smaller spectral extraction 
region, and the ability to pre-filter the particle background (see section 2; the Draco obser- 
vation was conducted in FAINT, rather than VFAINT, data-mode). We derive a 90% upper 
li mit at ~ 2.5 keV of 2.0 x 1 0"^ ph otons cm ^ s ^ in Draco, consistent with the results 



of iRiemer-Sorensen fc Hansen! ( 120091 ). Based on an estimated projected mass of ~ 6 x 10 



6 



'*http://cxc. harvard.edu/proposer/POG/ 
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Fig. 10. — (a) top: AC = 2.3, 4.61, and 9.21 confidence contours for tlie emission line 
parameters, from tlie Suzaku XISl spectrum of Ursa Minor, (b) bottom: XISl (red) 
and XIS03 (black) spectra and best fit models that include a ~ 2.4 keV line (the flux is 
2.5 X 10^^ photons cm~^ s"^ in the former, and in the latter). 
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Mq (IStrigari et al.ll2007al ). the expected flux is ~ 75% that in Willman 1 (Draco is twice as 



distant) - i.e. (2.5 ± 2) x 10 ^ photons cm ^ s 



4.2.2. Other Dark Matter Dominated Systems 

S ome of the strongest claimed hmits on sterile n eutrino parameters are derived from 
M31 (IWatson et al.l l2006l : iBoyarsky et al.l l2008l . 12010) ■ These depend on estimates of the 
dark matter content in regions where baryonic matter is significant, if not dominant. While 
there is a fairly strong consensus on the total mass in M31, there is a wide range of possible 
mass dec ompositions (primarily) into stella r disk, stellar bulge, and dark matter h alo com- 



ponents (IChemin. Carignan. fc Fosteiil2009t ICorbelli et al.ll2009l : ISaglia et al. 



2010) - wider 



than assumed in the papers referred to ab ove. The best-fits presented in ICorbelli et al. 
(120091 ) are obtained using the iBurkertI ( 1l995l ) dark matter density distribution function. The 
resulting projected dark matter mass of 2.9 x 10^ Mq (for models with constrained virial 
mass) in the 1.1-3 kpc annulus, demonstrates that a non-detection in M31 does not stand in 
contr adiction to a feature in Willman 1 with best-fit parameters we derive from the Chandra 
data (IKusenko and Loewensteinll2010l ). 



The diluting effects of X-ray emission associated with stellar and/or hot plasma presents 
a potential complication in systems with prominent baryonic components. For galaxies, 
we estimate the contribution from an old stellar population of mass M^, to the luminosity 
e mitted in an energy interval AE at 2.5 keV, using the power-law indices and normalizations 



m 



Revnivtsev et al.l (120081 ) . as 



LAB(2.5keV)/M, = 1.0 x 10^^[1 + 40(1 - /,es)] 



20 



(5) 



where the ratio AE/E is normalized to a typical CCD spectral resolution (at 2.5 keV) of 20. 
The first term accounts f or cataclysrnic var i able and coronally active binary stars, as seen in 
the Milky Way and M31 JlI fc WanglboO?! : Isogdan fc Gilfanovl 120081 ). The second accounts 
for low mass X-ray binaries, with /res the fraction of the LMXB flux that is resolved out 
(which depends on the angular resolution and exposure depth of the observation, as well 
as the LMXB luminosity function). The ratio of this term to the sterile neutrino radiative 
decay line emission from the dark matter mass Mdark is 



[0.111 + 4.44(1 - /,,J]^^r:^,— ^ 

Ma 



(6) 



E J W dark 

This emission, that may be augmented by discrete sources and hot gas associated with a 
more recent star formation, may be non-negligible where dark matter is subdominant or 
marginally dominant. 
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In galaxy clusters and groups, and massive elliptical galaxies, that are filled with hot 
(~ 1 — 10 keV) plasma enriched with the products of star formation, emission of a 2.5 keV 
sterile neutrino decay line may be overwhelmed by Ka emission lines at ~ 2.45 and ~ 2.6 
keV from He-like and H-like S, respectively. For example, in M87 the observed S XV line 
flux (IMatsumoto et al.lll996r) within 10' is ~ lOx that expected from the amount of dark 
matter ( Churazov et al.T2008 ) if it were composed of sterile neutrinos decaying at the rate 
inferred in Willman 1. 

Of course, the gravitational potential wells in dwarf spheroidals are far too shallow (their 
observed velocity dispersions corresponding to ~ 10 eV temperatures) to bind gas hot enough 
to coUisionally ionize S to this degree, and there are no significant sources of photoionization. 
If hot plasma were somehow responsible for the Willman 1 2.5 keV emission line, one would 
expect much stronger accompanying emission in Si-K and Fe-L shell emission that is not 
seen - unless the elemental abundances are highly and uniquely peculiar. When we attempt 
to fit the Willman 1 spectrum with a thermal plasma model (best fit temperature, kT ~ 2 
keV), we find that the 2.5 keV line may be accounted for only if the S/Fe abundance ratio 
exceeds 50 x the solar ratio. 



4-2.3. Other Constraints 



The limits based on Lyman-a forest data (IBoyarsky et al.ll2009l ) allow up to 40% of dark 
matter in the form of 5-keV sterile neutrinos produced via Dodelson-Widrow (DW) scenario 
(IDodelson fc Widrowl Il994j ). However, if there is any contribution from the Higgs boson 
decays, the resulting population of relic ste rile neutrinos has about three times smaller free- 
strea r ning length than th e DW population (jKusenkdl2006l : |Petrakill2008l : iPetraki fc Kusenko 
20081 : iBoyanovskyI l2008bl ) . The inclusion of the additional "cold" population of sterile neu- 
trinos weakens or eliminates the Ly-a bounds. Although a detailed numerical analysis of 
such mixed dark matter has yet to be performed, the approximate analytical r esults indicate 



that 5 keV sterile n eutrinos can account for 100% of dark matter in this case (jPetrakil 12008 
Bovanovskvll2008al lbll 



4- 2.4- Instrumental Features in Proximity 

Improper subtraction of the PB could conceivably produce artefacts in the subtracted 
spectrum that may be misinterpreted as a feature. We address this issue, as follows, by ap- 
plying an analysis procedure to the Willman 1 spectrum along the lines that we implemented 
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in analyzing the Suzaku dwarf spheroidal spectra (see Section 4.2.1, LKB). We consider the 
full unsubtracted, and unbinned, spectrum in the 2.2-5 keV band. The reduced-background 
PB spectrum (Section 2.4) is fit by a model consisting of a power-law and two Gaussian emis- 



sion l ines with energies corresponding to the brightest Au fluorescent features (jChartas et al. 



2000l ). and line- widths that are free to vary. The addition of a line at 2.51 keV does not 



improve the fit, and the AC = 2.7 line flux upper limit is 6 x 10~^ photons cm~^ s~^. 
If we then fit the total (unsubtracted) spectrum with the same baseline model, the ad- 
dition of a line at 2.51 keV improves the fit by AC = 5.8 and the best-fit line flux of 
3.5^2 5 X lO^^photons cm~^ s^^ (AC = 2.7 errors) is consistent with the estimate based on 
the PB-subtracted spectrum. If we allow the line energy to float, we find energy 2.57 ± 0.06 
keV and flux 4.2 ± 2.5 x 10~® photons cm~^ s~^ (AC = 4.6 errors), again consistent with 
the initial estimate. The improvement in fit over the baseline model is AC = 10.6. If we 
then go back and refit the PB spectrum with an additional line with energy fixed at 2.57 
keV, we find that the AC = 2.7 line flux upper limit is 1.6 x 10~^ photons cm~^ s~^, and 
the inclusion of the line does not improve the fit (AC = 1.6). These considerations indicate 
that the 2.5 keV feature is not an artefact of inaccurate PB subtraction. The unbinned PB 
and total spectra and best fit models are shown in Figure llab; the same comparison with 
the data rebinned (in the plot only) is shown along with the two best-fit models in Figure 
12. Confidence contours are plotted in Figure 13. 



4-2.5. Consistency, but not Corroboration 



Our re-examination of the other X-ray dwarf spheroidal datasets yields an intriguingly 
similar weak 2.5 keV feature in the Suzaku XISl spectrum of Ursa Minor, but fails to confirm 
the presence of an emission hne at the expected stre ngth in Suzaku and C handra data on 
Draco and the Suzaku XIS03 spectrum of Ursa Minor. iBoyarsky et al.l ( 120101 ) do not detect a 
feature in either the Sculptor or Fornax dwarf s pheroidals; altho ugh, it is notable that these 
systems have relatively low mass-to-light ratios (IWolf et al.ll2009[ ) - increasing the possibility 
of the presence of the diluting effects of an unresolved stellar component (the former is known 
to have a population of X-ray binaries; Maccarone et al. 2005) - a factor that is also an 
issue in M31 (see above). Detection of a 2.5 keV line from sterile neutrinos is impeded in the 
most massive dark-matter dominated systems, which are well known to retain hot metal- 
enriched (interstellar, intragroup, or intracluster) gas with line-emission from highly ionized 
S. Such sources of confusion are absent in the extreme baryon-poor faint and ultra-faint 
dwarf spheroidals, making them the cleanest targets for confirming or refuting the evidence 
for 5 keV sterile neutrino dark matter uncovered in Willman 1. 



- 24 - 



> 

ID 



0.04 



0.02 



5x10-3 



-5x10-- 




Energy (keV) 




Energy (keV) 

Fig. 11. — (a) top: Unbinned PB spectrum and best fit model (see text), (b) bottom: 
Same for total Willman 1 spectrum. Inclusion of an additional emission line, with energy 
and flux consistent with that obtained from the PB-subtracted data, significantly improves 
the fit in the total, but not the PB, spectrum. 
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That is, the current absence of other positive detections may, in every case, be accom- 
modated with confirmation of the Willman 1 provisional detection reported here - given 
the uncertainties in mass profiles and various other factors discussed in this section. Our 
measurement lies near the limit of faint source line fluxes attainable with the best current 
detectors and, even in the most promising targets, the possible presence of sterile neutrinos 
with the most suitable parameters in terms of viability as a dark matter candidate lie tanta- 
lizingly close to the edge of detectability. It is notable that the Willman 1 X-ray observation 
is unique in that it represents the only deep Chandra exposure of a very high mass-to-light 
system, so that all sources of competing 2.5 keV line and continuum emission of astrophysical 
origin are minimized. On the other hand, the dynamical nature of Willman 1 remains un- 
resolved and continues to be investigated (iGutowski et al.ll2010l : IWillmanll2010[ ). Additional 
observations of this - and/or other similar - dwarf spheroidals may be required to settle the 
question of whether the Willman 1 feature is real or spurious. 



5. Implications of a Detection 



One can take two different approaches to interpreting the spectral line from a decay- 
ing relic particle. If one assumes dark matter is solely composed of sterile neutrinos, then 
the number density of particles is determined by the particle mass and the mixing angle is 
inferred from the line flux. Alternativ ely, one can relax th e assumption equating the relic 
density and total dark matter density (IHinshaw et al.l 120091 ) - e.g., there may be additional 
contributing dark matter components. Regardless of any additional phy sics, sterile neutrinos 
of a g iven mass and mixing angle are produced by neutrino oscillations (IDodelson fc Widrow 
1994 ) at a rate that is calculable and cannot be reduced by any unknown high-scale physicsj^ 
The contours for the allowed ranges of mass and mixing angle are somewhat different de- 
pending on whether the relic sterile neutrino densit y is assumed to be equal to the measured 
cosmological density of dark matter (jKusenkdl2006l ) . We show the results of both approaches 
in Figure dU 

Interpreted as a sterile neutrino radiative decay, the emission line energy and 90% flux 
limits imply mst = 5 ± 0.2 keV, r^,^^^„ = (1.4 ± 1.1) x lO-'^'^{UM7/0.2)-\ and sin^ ^ = 
(7.8±6.1) X 10~^°(/stM7/0.2)~^, where the projected mass in the Chandraheam is normalized 
to the fiducial value of 2 x 10*^ M©. sin^ ^ = 1.7 x 10-^^{My/0.2)-^ represents an absolute 
lower limit to the mixing angle, regardless of the sterile neutrino production mechanism. If 



^The only caveat to this argu ment is the possib ihty of a non-standard cosmology with a low reheat 
temperature, as pointed out bv iGelmini et al.l (|2004l ). 
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Fig. 14. — For sterile neutrinos that make up 100% of dark matter, the allowed range for 
the mass and the mixing angle inferred from the Chandra data is shown by the solid con- 
tour. This range is consistent with the region corresponding to production by non-resonant 
oscillations alone that is circumscribed by the solid lines marked "dark matter (min)" and 
"dark matter (max)". If the primordial abundance of sterile neutrinos is not assumed, the 
allowed range for the mass and the mixing angle is shown by dotted, dashed, and dot-dashed 
contours corresponding to maximum, average, and minimum efficiency prod uction by neu 
trino o scillations. A projected mass within 55 pc of 2 x 10^ Mq is assumed (jStrigari et al. 
2008bl ). The inferred parameters are consistent with the pulsar kick mechanisms based on 
either resonant (MSW) or non- resonant production in the cooling neutron star (see text, 
Kusenko 2009). 
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there was some DW production, at least 9.5{Mj/0.2)~^% of the dark matter must have been 
produced in this manner (based on the minimum production and lower line flux limit). If 
the DW mechanism is responsible for all of the sterile neutrino production, sin^ 9 = (2.5 — 
15.7) X 10~^°(M7/0.2)~^/^, where the lower (upper) limit corresponds to the line flux lower 
(upper) limit and the maximum (minimum) production; and, fst > 0.31(M7/0.2)~^/^. The 
data is consistent with DW production and fst = 1 (Figure 9) for My = 0.019 — 0.77, where 
the lower (upper) limit is derived from the minimum (maximum) flux and production. The 
corresponding mixing angle for minimum (maximum) production is sin^ 9 = 17.7(3.6) x 10^^°. 
In fact our best fit fiux matches that expected for average DW production and fst = 1 for 
Ml = 0.18 (sin^ 9 = 8.2 x 10^^*^), which is close to our fiducial mass (Figure 9). 

The allowed region may be further circumscribed if combined with previously derived 
constraints. Limits have been derived using a variety of datasets (mostly based on Chandra 
or XMM-Newton data), employing a variety of data reduction, background mitigation, and 
statistical techniques to derive upper limits, taking diverse approaches (if any) to account for 
emission associated with baryons (see above), and depend on different methods for estimates 
(and associated uncertainties) of the dark matter content in the field-of-view. The limits 
based on Chandra deep field measure ments of the unresolve d X-ray background originating 
in the Milky way dark matter halo (lAbazajian et al.l 120071 ) are perhaps the most directly 
comparable in that deep Chandra ACIS-I data (some of which was taken in VFAINT mode) 
was utilized. These constraints are consistent with dark matter being composed of nist = 5 
keV sterile neutrinos produced by the DW mechanism. XMM-Newton observations of diffuse 
emission in the M31 halo are formally inconsistent with our best-fit fiux under the same 
assumptions. However, taking the uncertainties into account, an ample region of concordance 
is retained that includes parameters where fst = 1 and non-resonant oscillations are the 
dominant, or sole, mechanism of sterile neutrino production (Section 4.2). 

This interpretation has immediate implications for particle physics and astrophysics. In 
particular, supernova physics is altered dramatically by the emission of a 5 keV sterile neu- 
trino that has a mixing angle sin 9 ~ 10~^ — 10~^. First, the emission is anisot ropic, with the 
direc t ion of anisotropy s et by the electron spins polarized in the magnetic field ( iKusenko fc Segre 
19971 : iFuUer et al.l |2003| ) . While the overall energy removed by sterile neutrinos with 5 ke V 
mass and sin 9 ~ 10~^ — 10~^ is small, the anisotropy in the emission of sterile neutrinos 
plays a central role in generating the supernova asymmetries. Second, the lepton number, 
energy, and entropy transport in the cooling neutron star are affected, and an increase in the 
electron neutrino luminosity can pr eheat the material in f ront of the shock wave, adding to 
the energy of the overall explosion (IHidaka fc FuUen 120071 ). Third, the neutrino-driven kick 
enhances the co nvection in front of the neutron star, which causes some energy increase of 
the shock wave (IFryer fc Kusenkdl2006l ). The numerical simulations of this effect show the 
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formatio n of asymmetric 



motion (IFryer fc Kusenko 



ets w ith the stronger jet pointing in the direction of the pulsar 



2006f ) 



From the point of view of particle physics, the 5 keV sterile neutrino with a small mixing 
angle can be accommodat ed in a seesaw m ass matrix that satisfies all the experimental and 
observational constraints (IKusenkd 120091 ) . The values inferred from the Chandra observa- 
tions allow for all dark matter to be produced via neutrino oscillations, so that no additional 
production mechanism is necessary. This, however, does not preclude the 5 keV mass from 
being generated via the Higgs mechanism by a vacuum expectation value of an electroweak- 
scale singlet Higg s boson, which may be discoverable by up coming experiments at the Large 
Hadron Collider kusenkol I2OO6I : IPetraki fc Kusenkol boosh . 



The X-rays emitted by dark matter decay during the "dark ages" could produce a strong 
enough ionization to enhance the formati on of molecular hydrogen, which, in turn, could 
speed up the formation of the first stars (IBiermann fc Kusenkol l2006l ). The corresponding 
predictions for 21 cm and other observa tions depend on the knowledge of the particle mass 
and mixing angle (jStasielak et al.l 120071 ) that are greatly refined if the dark matter particle 
has the mass and mixing inferred from Chandra observations of Willman 1. 

Finally, there is no doubt that observational cosmology will take advantage of our finding 
if it is confirmed. If dark matter emits a narrow spectral line in X-rays, one can hope to map 
out the three-dimensional dark matter distribution, including the redshift, using future X- 
ray telescopes such as and Such data could be used to study the expansion 

of the universe on cosmological distance scales. 



6. Summary and Concluding Remarks 

We presented results of a search for a radiatively decaying dark matter emission line 
in the Chandra spectrum of the ultra-faint dwarf spheroidal galaxy Willman 1. Ax^ = 9.2 
upper limits on the line flux over the 0.4-7 keV bandpass were translated to an allowed region 
in the sterile neutrino mass-mixing angle plane for the case where all of the dark matter is 
composed of sterile neutrinos regardless of how they are produced, and in the more general 
case where their abundance is not assumed but is calculated under the assumptions that 
minimize it. These constraints are consistent with and therefore, given the uncertainties in 
mass determination, reinforce those derived from previous X-ray analysis of dwarf spheroidal 



^http://heasarc. gsfc.nasa.gov/docs/astroh/, 
^http:/ /ixo. gsfc.nasa.gov/ 



http://astro-h.isas.jaxa.jp/ 
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galaxies (LKB, Riemer-Sorensen & Hansen 2009). 

We uncovered evidence of an emission line in the Willman 1 X-ray spectrum with energy 
2.51 ± 0.07(0.11) keV and flux [3.53 ± 1.95(2.77)] x 10"^ photons cm'^ s-^ where the errors 
are 68% (90%) confidence limits. This was based on an excess to the continuum measured by 
fitting the particle-background-subtracted source spectrum, where the VFAINT data-mode 
quiescent background reduction technique was applied. Consistent parameters were found 
using maximum likelihood fitting to the unbinned, unsubtracted spectrum. Re-examination 
of Suzaku data of the Ursa Minor dwarf spheroidal, and Chandra data on the Draco dwarf 
spheroidal, reveal a possible confirmation in the Suzaku XlSl spectrum of Ursa Minor, but 
not in any of the other spectra. Evidence of dark matter radiative decay in Willman 1 
remains provisional due to the low level of significance, but uncontradicted. In the short 
term, additional X-ray observations of SDSS dwarf spheroidals are needed, while future high 
energy resolution X-ray spectroscopy should prove definitive. 

We derive the range of allowed sterile neutrino mass and mixing angle implied by the 
90% confidence hmits of line energy and flux, assuming that either (1) dark matter solely 
consists of sterile neutrinos regardless of how they are produced, or that (2) they are solely 
produced by neutrino oscillations regardless of whether the abundance matches that of dark 
matter as inferred from cosmology. Remarkably, the two allowed regions strongly overlap 
- that is, the observations are consistent with the narrow region in the mass-mixing-angle 
plane where neutrino oscillations may produce all of the dark matter (although it allows 
for some non-DW production). Moreover, the inferred parameters are consistent with the 
pulsar kick mechanisms based on either resonant or non-resonant sterile neutrino production 
in cooling neutron stars. This unlikely concordance intensifies the importance and urgency 
of additional investigation. 
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